Determination of time-varying pressure field from phase contrast MRI data

Methods We propose a numerical framework for determining pressure field in large arteries, the main components of which are: axi-symmetric 1-D unsteady wave-propagation model with elastic walls and an optimization framework for model personalization via parameter estimation from measured flow and anatomical data (Fig 1). 4D PC MRI data was acquired using a multi-slice phase contrast FLASH cine sequence with 3d velocity encoding (Siemens Magnetom Verio, temporal resolution 52 ms, matrix 100x192, FOV 319x219 mm, slice thickness 5 mm) in a sagittal to coronal slice orientation. Siemens 4D Flow WIP software was used for background phase, motion and venc anti-aliasing corrections, extracting vessel models including centerline trees and lumen boundary, and for automated extraction of flow profiles and cross-sectional areas. We use a linear elastic model for the vessel wall, with an axially varying stiffness. The parameters of the material model are estimated by a numerical optimization procedure, wherein the simulated deformation is fitted to the observed deformation of the vessel wall extracted from the MR data. Unsteady flow rate at the aortic inlet is extracted from PC-MRI data and used as the inlet boundary condition. To obtain physiologically valid mean and pulse pressure values and the wave propagation effects, we use a structured tree outflow boundary condition, wherein a morphological or fractal vessel tree is introduced to model the downstream behavior at the outlets of the vessel tree. The parameters for the outlet boundary condition are then independently perturbed at each individual outlet in order to match the measured unsteady outflow rate from PC-MRI data. This is achieved by adapting the final radius at which the asymmetric binary tree is terminated. The hyperbolic PDEs are solved using finite differences based on 2nd-order Lax-Wendroff method. Simulations converge after 3 cycles, with a maximum relative difference of 1% in pressure and flow between 2 consecutive cycles.